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(54) Procede de localisation spatio-temporelle, de radio-mobiles en milieu urbain 

(57) Le procede consiste : 


a calculer (1 ) la r6ponse impulsionnelle du canal de 
propagation entre le radio-mobile et la station emet- 
trtce/r^ceptrice, 

a estimer (2) k partir de la reponse impulsionnelle 
les couples (angles, retard) du trajet de propagation 
direct et des multitrajets r6sultats de reflexions sur 
des obstacles situ6s dans le chanp de propagation 
des ondes entre la station de base et le radio- 
mobile. 

a calculer (2) pour chaque trajet les attenuations Oj. 
et k calculer (3) la distance d s6parant la station 
emettrice/r6ceptrlce du radio-mobile en consid6- 
rant seulement les retards et les attenuations du 
trajet direct et d'un multitrajet quelconque. 

Applications : Localisation etsuivi de radio-mobiles. 
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Description 


[0001 1 La presente invention concerne un proc6d6 de localisation de radio-mobiles en milieu urbain. 
r00O2] Elle s'applique notamment k la localisation et au suivi de radio-mobiles dans un r6seau GSM pour antiaper 
5 leur passage dune cellule k la suivante et amfeliorer la reception des stations de base par la connaissance de la posi- 
tion des radio-mobiles 6metteurs. 

[0003] Dans des configurations simples de propagation, la localisation de radio-mobiles est traditionnellement obte- 
nue a I'aide de radio-goniom6tres dispos6s en des lieux determines rfun territoire. La position d'un radio-mobile peut 
alors 6tre determin6e par triangulation en reportant sur une carte geographique les lieux oil sort srtu6s les radiogonio- 

70 metres ainsi que les directions d'arriv6e sur ceux d des ondes radio6lectriques 6mises par le radio-mobile. 

100041 Cependant cette m6thode qui repose sur Texploitation des seuls chemins de propagation directs apparait dil- 
ficilement applicable dans le cas ou les ondes empruntent plusieurs trajets et est 6galemem lourde de mise en oeuvre 
car elle necesslte plusieurs stations de reception. Ce type de propagation est notamment celui qui est rencontrfe en 
milieu urbain dans les r6seaux de communication radio-mobiles cellulaires. ou le spectre des communications est situ6 

,5 dans une gamme d'ultra hautes frequences s'^tendant typiquement de 30 MHz k 3 GHz pour le r6seau GSM. Pour 
cette gamme de frequences la propagation des ondes radlo6lectriques entre un 6metteur et un r6cepteur peut 
emprelnter des trajets multiples dus aux reflexions et aux diffractions sur des obstacles interposes sur le chemin de pro- 
pagation direct entre I'emetteur et le r6cepteur. La gamme des frequences utilisees est aussi la cause de trajets diffus, 
aussi bien temporellement que spatialement, qui degradent les perlormances des liaisons. 

20 [0005] Le but de I'invention est de pallier les inconvenlents precites. ... 

[O0O6] A cet effet. I'invention a pour objet un proc6d6 de localisation spatio-temporelle d'un radio-mobile relativement 
k une station de base emettrice/r6ceptrice caract6rise en ce qu'il consiste : 

- k calculer la r6ponse impulsionnelle du canal de propagation entre le radio-mobile et la station emettrice/recep- 

25 trice 

. k e^imer k partir de la r6ponse irrpulsionnelle les couples (angles, retard) du trajet de propagation direct et des 
multitrajets resultats de reflexions sur des obstacles situ6s dans le champ de propagation des ondes ertre la sta- 
tion de base et le radio-mobile. 

- a calculer l attenuation Qj de I'onde propagee par chaque trajet. 

30 - et a calculer la distance d s6parart la station emettrice/receptrice du radio-mobile en considerant seulemert les 
retards et les attenuations du trajet direct et d'un trajet non direct quelconque. 

[0007] D-autres caracteristiques et avantages de Unvention apparaitront k Taide de la description qui suit faite en 
regard des dessins annexes qui repr6sentent : 

- la figure 1 le precede selbn I'invention sous la forme d'un organigriamme, 

- la figure 2. un modeie de reception monocapteur, 

- la figure 3, une illustration d'une reception multicapteur par une antenne circulaire. 

- les figures 4 et 5. une illustration d'une'transmission k deux bu trols trajets en milieu urbain. 

[00081 Le precede selon I'invention se deroule de la fagon representee k la figure 1 suivart trols etapes principales. 
Dans une premiere etape. representee en 1 & la figure 1 . il consiste k calculer la r6ponse impulsionnelle fi(t) du canal 
de propagation pour estlmer k l'6tape 2 les angles e„ et les retards t„, des diff6rents trajets ainsi que les attenuations 
ou des ondes propag6es sur chacun des trajets et enfin k calculer k I'etape 3. k partir des retards t„, calcuies et de la 
45 reponse Impulsionnelle h(t) du canal, la distance d s6parart un radio-mobile de la station d'6missionyr6cepton. 

[00091 Le calcul de la r6ponse impulsionnelle du canal qui est effectu6 k retape 1 a lieu sur le signal obtenu k la sortie 

d'une antenne multi-capteurs. . . ^ \. , ^- w 

[001 01 En presence de plusieurs trajets M d'6mission bande etroite de radio-mobiles, le signal obtenu k la sortie d un 
capteur n de l antenne muW-capteurs est egal k la somme des M signaux s„(t) 6mis dans les directions retard6s de 
50 T et multiplies par la reponse a,{e^ du capteur i dans les directions 6^ k taquelle s ajoute un bruit bj{t). Dans I'hypo- 
these d une transmission num6rique k bande etroHe de frequence, une representation de la sortie vectonelle de 
I'artenne compos6e de fM capteurs s'obtient en concatenant la sortie Xj(t) de chaque capteur, chaque sortie Xi(t) ven- 
f iant la relation : 
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(t)= £s„(t).ai(e (1) 


m=»1 
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[001 1 ] En introduisant un vecteur directionnel B{e^) tel que 

a(ej = [ai(ej.....aN(ej]^ (2) 
et en concatenant ta sortie Xj(t) de chaque capteur. on obtient un vecteur x(t) tel que : 

M 

X(t) = [X-,(t). Xj(t). X^O)]^ = £ 9(Qm)SmO)+b(t) (3) 

[0012] En rassemblant les vecteurs directionnels 3(6^) dans une nnatrice 

A = [a(ei),a(e2) sOn)] 

la relation precedente s'ecrit encore : 

x{t) = As(t)-hb(t) (4) 


avecs(t) = [s,(t).S2(t)- • 

ou Si est le signal du 1®^ trajet. 

S2 est le signal du 2'®"'® trajet. 

Sm est le signal du M'®"^® trajet. 

[001 3] Suivant le modele de propagation qui est represents schematiquement a la figure 2, te signal x(t) qui est 
obtenu a la sortie d un capteur de reception est un signal de la forme : 

x(t)= £a^.e'*ra(0J.s(t-Tj + b(t) (5) 

m=1 


OU designe I'attenuation du m'^"^® trajet 

est le retard de groupe du m'®"^® trajet 
M est le nombre de trajets 
et <|)rT, est le d6phasage du m'®"^ trajet. 

soit encore sous la forme d'une convolution du signal 6mis s(t) avec la r^ponse impulsionnelle h(t) du canal : 

soitx„(t) = h„(t)-s(t) + b„(t) (6) 


oil bn(t) reprSsente I'enveloppe complexe du bruit additif . 
[0014] Dans le cas d'une reception sur I'antenne composSe de multi-capteurs, comme cela est represents k la f igure 
3, la reponse impulsionnelle du canal est modif i6e par la reponse de I'antenne suivant les directions d'arriv6e. Cette 
rSponse est de la forme : 


hn(t)= X:^n(er..Ym)«me"'"6(t-xJ (7) 
m=1 
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en supposant te canal ideal, c'est-a-dire de bande infinie. 
[001 5] En mode 6chantillonne. la r§ponse du canal de bande B prend la forme : 


h nM = F ^ ^ a J e ^) a ^ e ' * Sin c((k -F^xJ B) 


(8) 


m=1 


avec k e Z. i'espace des nombres entiers. 


10 


sin(7i(k-FeT^)B) 
et sin c((k - F ,x^) B) = ^(k.F^,^)B 


(9) 


15 est la frequence d'6chantillonnage 

B est la bande en frequence r^duite du signal. 

[001 6] Cette r^ponse peut §tre estim^e par un f iltre causal, a reponse impuisionnelle f inie de longueur L et de coef- 
ficients tels que : 


20 


h n[k] = £ a e ^) a ^ e '* sin c((k - F ^) B) 


(10) 


25 


pour ke {0 ... L} 

et hjkl=0pourk (Oouk )L 


30 


Le vecteur t3-,correspondants*6crithn = lhn(0). • • - .h^CL)]^ 


(11) 


35 


[0017] Suivant un premier mode de realisation de Tinvention. I'estimation des angles 6^ et des retards de r^tape 
2. est executee de fagon disjointe. en utilisant la relation matricielle suivante : 

H = [hi.h2. .bNl^l^'^A^ (12) 


avec : 


40 


45 


M = 


sinc(-TiFeB) sinc(-X2FeB) 
sinc((1-TiFe)B) sinc((l- X2Fe)B) 


sinc(-TMFeB) 
sinc((1-TM'^e)B) 


sinc((L-tiFe)B) sinc((L- TzFg) B) sinc((L- T^/,Fe)B) 


(13) 
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a-i 0 

0 OL2 


10 


0 ••• 0 ttN 
A = [ai(0i)...aN(eN)] 


(14) 
(15) 
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[0018] L'estimation des parametres temporels s'obtient k partir de la relation (12) en estimant un projecteur 

sur I'espace temporel HH^ et en appliquant soil la nn6thode connue sous I'abreviation anglo-saxonne MUSIC, 
de < (Mutti Signal Classification >>en cherchant les maxima de la quantite 

Pt= r-^ (16) 

m(t)^P^^t!I?(x) 


OU 


m(T) = [sinc(-TFeB). .... sinc((L - TFe)B)]* 

est une colonne de la matrice M(13), soit la m^thode deterministe connue sous le nom de ( < maximunh de vraisem- 
blance >) en utilisant un crit^re de la forme 

J(ti, .... tm)= tr -M (MtM)-^IVlt) HH^) (17) 


40 [001 9] Une description de la m6thode MUSIC peut etre trouv6e dans I'article de M. Ralph O. Schmitt ayant pour titre 
< < Multiple source DF signal processing : An experimental system ) > publi6e dans IEEE transaction on antennas and 
propagation AP34 N° 3 March 1986. 

[0020] La determination des M couples (Grr,, t^) permet de localiser les radio-mobiles. Dans ie cas ou. comme cela 
est represente k la figure 4. seulement deux trajets, un trajet direct et un trajet comportant une reflexion sont selection- 
45 n6s. la determination des couples a lieu en consid6rant quatre parametres , . 02, t2 qui se ram^nent ^ la connais- 
sance de 0i , 62 et des retards , t2. En appelant d la distance du trajet direct qui s6pare I'^metteur radio-mobile 4 d'un 
r6cepteur 5 etd-j, d2 les distances du trajet indirect s6parant respectivemerrt le r6cepteur 5 d'un obstacle 6 sur lequel 
Tonde indirecte est refl6chie et Temettetur 4 de Tobstacle 6, la distance d est fournie par la relation : 

. (T2-T,)C 


55 ou et 02 designent les attenuations du premier et du deuxieme trajet et \i la perte due ^ la reflexion du second 

trajet Les attenuations et 02 sont obtenues k partir des couples (0^. ; Og. 13) et de la relation matricielle 
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a=M*H A** = 


ai 0 
0 a2j 


(19) 


ou correspond k I'attenuation du i^""® trajet 
[0021] Comme la perte de puissance dans une propagation libra se fait en 


10 


,2 


15 


ou r est la distance parcourue, et U2 sont li^s par la relation : 


(20) 


25 


_ (T2-T^)^c^+2d(Tg-T^)c 
avecdi - 2[(-c2-Xi)c+d-dcos(e2-e^)] 


etdg = j6^+6^^'26di^cos(Q2 ' ©i) 


(21) 


(22) 
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[0022] S*il y a plus de deux trajets, la figure 5 montre que Tinformation apportee par tout trajet suppl6mentaire se 
resume a deux parametres. et les nouvelles inconnues sont ^alement au nombre de deux. En consequence, une 
detection de plus de deux trajets n'a pas d'inter§t dans le sens ou elle ne supprime pas les inconnues. II suffit dans ces 
conditions de limiter les traitements en considerant seulement le trajet direct et un trajet autre quelconque. Connaissant 
les angles d'arrivee des difierents trajets, il est possible, gr§ce a une carte topographique d^taillee en trois dimensions 
de trouver les points de reflexion les plus probables et ainsi de remonter a la source d'6mission. 
[0023] D'autres modes d'estimation des parametres sont 6galement possibles. 

[0024] Lestimation disjointe des parametres d^crite pr6cedemment a I'avantage de r^duire Tespace de recherche. 
Cependant, dans cette methode . les angles et les retards ne sont pas associ6s si bien qu'elle souffre d'un manque de 
resolution pour estimer les retards proches. 

[0025] D'autres modes de mise en oeuvre de Tinvention permettent d'associer les parametres en estimant la matrice 
des attenuations ii avec les diff^rentes combinaisons (e^. t^) possibles en cherchant les couples rendant la matrice 
a la plus diagonale possible. 

[0026] Selon un deuxieme mode de mise en oeuvre de I'invention, le proc6d6 op6re par estimations successives en 
estimant successivement les angles e et les retards x tout en les associant. Pour cela. deux m6thodes peuvent etre uti- 
lisees. La premiere consiste a estimer les parametres temporels t^, connaissant les angles 6^. 

[0027] En estimant des angles 6 ^ paria m6thode MUSIC ou la methode du maximum de vraisemblance I'estimation 
des retards est obtenue en estimant la matrice M par la relation matricielle 


M = HA» 


(23) 


50 


[0028] Chaque colonne de la matrice M forme un projecteur permettant de trouver le retard correspondant en 
appliquant I'algorrthme MUSIC pour une seule source d'6mission. 
55 [0029] Le critdre utilis6 est de la forme : 
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P(x) 


_ m(T)^P^. m(T) 

a}^(T)!i}(T) 

!D(ii)m^(ti) 


(26) 


10 


[D^(tj)m(Tj) 
avec m(x) tel que M = [m(ti) ... rrKx^^)] 
sinc(-T Fe B) 


et rn(T) = 


sinc((L-TFe)B) 


[0030] L'estimation des retards peut 6galement dtre obtenue par correlation. En consid6rant que le signal re$u sur 
un capteur de M sources est de la forme : 


20 


M 


X(t)= 52^(®m)Sm(t) + b(t) 


(25) 


25 


30 


et connaissant les angles 6^. il est possible de considerer un signal : 

s[k] = A*^x[k] 

qui est form6 de M versions du signal emis retardees tel que : 


35 


40 



■s(k-TiFe)" 

s[k] = 

s(k-X2Fe) 


s(k-TM Fe) 


(26) 


[0031] L'autocorr6lation d'une des composantes de s [k] permet d'obtenir une r6f6rence temporelle et I'lntercorr^la- 
tion de cette composante avec les autres composantes donne une estimation des retards relatifs. 
45 [0032] De fagon similaire k la m6thode pr6c6dente, il est aussi possible d*estimer les param6tres spatiaux 0^^ con- 
naissant les retards x^. Ceux-ci sont obtenus en calculant equation matridelle : 

=(m^m]"^M'^H (27) 


[0033] Chacune des m6thodes utilis6es dans ce deuxi^me mode de realisation du proc6d6 selon Tinvention permet 
apr6s avoir estim§ les angles ou les retards de trouver les autres inconnues par M recherches en une dimension. Elles 
55 ont pour avantage de repousser certaines ambiguTt6s. mais la seconde recherche apparaTt tr6s d^pendante de la qua- 
lite de la premiere estimation. Une alternative serait d'effectuer, comme cela a et6 decr'rt par M. J. GROUFFAUD dans 
son memoire de these de Doctoral de TEcole Normale Superieure de Cachan de juin 1997 intituiee < ( Identification 
spatio-tennporelle de canaux de propagation k trajets multiples ) > , des recherches successives en alternant les m6tho- 
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des pour affiner ta precision de restimation des parametres. 

[00341 Enfin, suivant un troisieme mode de realisation du proc^de selon Tinvention. I'estimation des parametres 0^ 
et a lieu de fagon conjointe en appliquant des m^thodes de determination des r6ponses impulsionnelles fondles sur 
une"d6composition en sous-espaces, un sous-espace signal de dimension M 6gale le nombrede sources et un sous- 
espace bruit. conpl§mentaire du sous-espace signal de dimension N(L + 1)-M. ou en appliquant la methode du maxi- 
mum de vraisemWance. 

[0035] Pour estimer conjointement les angles et les retards par la m6thode des sous-espaces. un projecteur Prh est 
cr66 a partir de la matrice d'autocorr^lation des filtres 


15 


20 


25 


[0036] L'estimation h du vecteur h = [h JO], h J1] h,[L] hN[0]. hM[1] h^[L]]' a lieu en recherchant les 

maxima de la quantite 


ou g(e.T) = a(e)«>m(T) el 


M 


m=1 


avec 


30 n - rr o*4*m 


[0037] Dans la recherche des parametres par la methode du maximum de vraisemblance les parametres cherches 
35 sont les minimum du critere. 

ou G = [g(ei.ti)...g(em.Tm)] 


45 [0038] Cette recherche a lieu dans Tespace de dimension M x 2 des parametres {Q^, t^. • • • , e^, x^). 

[0039] L'optimisation des difterents modes de realisation du procede selon Tinvention peut avoir lieu en appliquant 
une methode connue d'optimisation du type Gauss-Nevyrton. Apres avoir trouve un point d'initialisation par les m6thodes 
decrites precedentes, une optimisation des angles 6^ peut avoir lieu en deroulant I'algorithme suivant : 

50 Initialisation : 0© ®init 

Iteration k : 

Calcul du gradient v^-i et du hessien H^.i au point 6 
Determiner 

ek = ek-i-Hj^liVK_i 


8 


EP0 926 510 A1 


Test de la convergence 

si D < £ STOP sinon k = k + 1 

[0040] Le meme type d'optimisation peut etre utilise pour optimiser la recherche des retards Tj^. 
[0041] Un autre nnode d'optimisation peut etre mis en oeuvre par une version approch6e du maximum de vraisem- 
blable suivant par exemple la m§thode d^crite dans IEEE transaction on signal processing, VoL 39, N11, Novembre 
1 991 des auteurs M. Viberg. B. Ottersten. T. Kailath ayant pour titre : ( < Detection and estimation in sensor arrays using 
weighted subspace fitting ) >. 
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ANNEXE A 
GLOSSAIRE 

A.1. Conventions standards : 

X* : conjuguee de x 

x^ : transposee du vecteur x 

x^ : transposee conjuguee du vecteur x 

X* : transposee de la matrice X 

X^ : transposee conjuguee de la matrice X 

X*^ : inverse de X 

X** : pseudo-inverse de X : X* = {x^X^ X^ 

tr(X) : trace de X 

E[.] : esperance mathematique 

X, X : valeurs estimees 

X* y : convolution des signaux x et y 

5(t) : distribution de Dirac 

sinc() : fonction sinuscardinal 

® : produit de Kronecker 

A.2. Constantes entieres : 

N :nombre de capteurs composant Tantenne 

M : nombre de trajets 

L : longueur de la reponse impulsionnelle 

K : nombre d'echantillons temporels pour composer robservation 
spatio-temporelle. 

A.3. Variables scalaires : 

Xn(-) :observation sur le n'^™capleur 
s() : signal emis par la source 
b() : bruit 

: variance du bruit 
hnO : reponse impulsionnelle du n'^""* canal 
F, : frequence d'echantillonnage 
Tb : periode d'echantillonnage 
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am : amplitude du m'*™ trajet 
(t)m - dephasage du m**'"* trajet 
Qm : azimut du m'^""* trajet 

: site du m'*™ trajet 
Tm : retard de groupe du m'^™ trajet 
an(0) : gain du n'^™capteur dans la direction (9) 
J() : critere a optimiser 

A.4. Variables vectorielles et matricielles : 

X [Observation spatiale sur I'antenne. de longueur N 

y :observation spatio-temporelle sur I'antenne, de longueur NK 

h : vecteur des reponses impulsionnelles, de longueur N(L+1) 

a(9) : vecteur directionnel spatial, de longueur N 

rn(T) ; vecteur tempore!, de longueur L+1 

Iq : matrice identite de taille qxq 

Rx : matrice de covariance du vecteur x 

Ry : matrice de covariance du vecteur y 

Rh : matrice de covariance du vecteur h 

M : matrice des M vecteurs temporels : M = [m(Ti), mC'^M)] 

A : matrice des M vecteurs directionnels : A = [a(0-|), a(9j^)] 

fl : matrice des attenuations des M trajets 

H : matrice des filtres mis en colonnes. de taille (L+1)xN 

G : matrice formee des reponses impulsionnelles 


EP0 926 510 A1 

ANNEXE B 

LE MAXIMUM DE VRAISEMBLANCE DETERMINISTE 
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SO 


En considerant un signal de la forme : 
Yk = A(G)Sk +b|, 

et le bruit comme elant gaussien. centre de variance a^, la 
distribution de probabilite du bruit est de la forme : 

f(bk)--r7-:rTexp ° 


La distribution de probabilite des donnees qui est a maximiser est 
donnee par la relation ; 


1 


la relation 


La vraisemblance des donnees pour le vecteur des donnees 
concatenees est alors : 

|-K'<2 Xk)=nn^^^P 

k=in'^|a^l 

Pour maximiser la vraisemblance cela est equivalent a minimiser 

K 

-iog(L(xi, X2 Xk)) = -MKIoga^ + Xl^'k " ^Skl 

k=1 

La minimisation de ce critere par rapport a 9 amene a minimiser : 

EK = Zl'^k -ASkI 
k=1 

^Sk=(A^A) A^k 
d'ou : 


K 

k=ll 

comme 


-aJa+a) 


K 


k=1 


KP="k"k = tr(uJuk)=tr(ukuj) 


55 
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devient ; 



Revendications 

1. Precede de localisation spatio-temporelle d'un radio-mobile relativement k une station de base 6mettrice/r6cep- 
trice caract6rls6 en ce qu'il consiste : 

a calculer (1) la r6ponse impulsionnelle du canal de propagation entre le radio-mobile et la station emet- 
trice/receptrice. 

a estimer (2) a partir de la r6ponse impulsionnel!e les couples (angles, retard) du trajet de propagation direct 
et des multitrajets r6sultats de reflexions sur des obstacles situes dans le champ de propagation des ondes 
entre la station de base et le radio-mobile, 
a calculer (2) pour chaque trajet les attenuations <x^, 

et k calculer (3) la distance d s^parant la station 6mettrice/receptrice du radio-mobile en considerant seule- 
ment les retards et les attenuations du trajet direct et d un multitrajet quelconque. 

2. Procede selon la revendication 1 caracterise en ce qu'il consiste a determiner la distance d par la relation : 


(■t2-'^l)C 



ou T2 - est I'ecart des temps de propagation entre le trajet direct et le trajet reflechi, 
c est la ceierit6 des ondes transmises sur le canal de propagation, 


2:1 
a2 


est le rapport des attenuations des propagations sur le trajet direct et le second trajet, 
et \i ta partie due k la reflexion de Tonde sur le second trajet. 


3. Procede selon la revendication 2 caracterise en ce qu'il consiste pour calculer (2) les att6nuations et des 
ondes sur le trajet a calculer une matrice diagonale Q det inie par le produit matriciel : 


ou 



est la matrice des coefficients de la reponse inpulsionnelle de canal. 
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H = [hi.b2.".hN] 

et A ' est la matrice transposfee des vecteurs directionnels a„ (e„, y„) correspondant aux multtrajets. 

Proc6d6 selon Pune quelconque des revendications 1 k 3. caract6ris6 en ce qu'il consiste a determiner la reponse 
impulsionnelle du canal (1) au moyen d'une antenne multi-capteurs coup6 k un filtre causal, k reponse impulsion- 
nelle f Inie. 

Proced6 selon l une quelconque des revendications 1 a 4. caract6ris6 en ce qu'il consiste k effectuer une estima- 
tion des angles etdes retards en estimant les paramfetres temporels k partir de la connaissance des angles 6^. 

Proc6de selon l une quelconque des revendications 1 k 4. caract6ris6 en ce qu'il consiste k effectuer une estima- 
tion des angles et des retards en estimant les paramStres spatiaux 6^ k partir des retards 

Proc6de selon l une quelconque des revendications 1 k 4, caract6rise en ce qu'il consiste k effectuer une estima- 
tion conjointe des angles On, et des retards Xm- 

Proced6 selon Tune quelconque des revendications 1 a 7. caract6ris6 en ce qu'il consiste k estimer les 
parametres temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de ralgorithme connu sous 
le nom de (( MUSIC )>. 

Proc6d6 selon I'une quelconque des revendications 1 a 7. caract6ris6 en ce qu'il consiste a estimer les parametres 
temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de I'algorithme d6terministe connu sous 
la denomination de ( < maximum de vraisemblance ) >. 
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